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Abstract 

In this paper we investigate the one-dimensional parabolic-parabolic Patlak-Keller-Segel model 
of chemotaxis. For the case when the diffusion coefficient of chemical substance is equal to two, in 
terms of travelling wave variables the reduced system appears integrable and allows the analytical 
solution. We obtain the exact soliton solutions, one of which is exactly the one-soliton solution of 
the Korteweg-de Vries equation. 
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I. INTRODUCTION 


In this article we consider the mathematical model of chemotaxis introdnced by Patlak in 
the 1950s [l| and by Keller and Segel in the 1970s jsl-Q]- Chemotaxis, or the directed move¬ 
ment of cells (bacteria or other organisms) along the gradient of concentration of chemical 
irritants, attracts a great interest. At present there is a nnmber of very interesting results 
concerning the existence and properties of both regular and blow-up solutions of the above 
system in different dimensions The simplihed model has the form: 

Ut — V{SiVu — r]iuVv) = 0 ^ ^ 

vt - - r]2U = 0 , 

where u{t,x) is the cells density, v{t,x) is the concentration of the chemical substance. The 
positive constants hi, § 2 , hi) h 2 ai’e the cells and chemical substance diffusion coefficients, 
the chemotaxis constant and the rate of cells chemical production respectively. After the 
replacement f —)■ hit and v the coefficients in the first equation in (1) become equal 

to unity. 


II. EXACT SOLUTION 


( 2 ) 


We investigate the following one-dimensional parabolic-parabolic model; 

Ut Uxx T (llUx)® 0 

Vt - aVxx - I3u = 0, 

where x G 3fJ,f > 0, n = u{x,t), v = v{x,t), a,/? are positive constants. This system 
is invariant under the one-parametric group of translation (x, t, u, n) —>• (x -|- ce, t + e,u,v). 
Rewriting (2) in terms of travelling wave variables y = x—ctwe obtain the system of ordinary 
differential equations (ODEs); one can see, that the first equation may be integrated, and 
on can put = 1. The reduced system has the form: 

Uy + CU — UVy -|- A = 0 

( 3 ) 

aVyy + CVy -f M = 0, 

where u = uiy}, v = v{y), A = const. Using the Kovalevskaya-Gambier method (see j^), 
improved in [10], [ll|, we can see, that this system possess the Painleve property only if 
a = 2. Let us focus on this case. 
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It is convenient to solve Eqs.(3) in terms of variable 


K, cy 

c 


( 4 ) 


where k > 0 is an arbitrary constant. Then for v and u we obtain the solutions in the form: 


z; = — In 


K 


zZ^iz) 


( 5 ) 


U = ( 1 — 


'2 ' - where z/2 = 1 - A 


The function Z^{z) satishes the modified Bessel’s equation and can be present as a linear 
combination of Infeld’s and Macdonald’s functions. 

Using the series expansion of the Bessel’s functions, as well as theirs asymptotic behaviour 


12|, one may obtain the following asymptotic forms for and Uy{z)\ 


z —)■ oo : —)■ 0; Uy{z) —)■ 0. 


( 6 ) 


z -)■ 0 : 




oo, 


K 


Stt 


0 < < 


2 ’ 


IcIU^ (7r + 2)2’ 

0 , 


= 


2 ’ 


1/ > 


2 ’ 


( 7 ) 


u,y z) c ly 


( 8 ) 


where the expression for ^ — 2 with the Eqs.(9) below. 
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Consider now the class of solutions with half-integer index ly = n-h when Z^(z) can be 
expressed in hyperbolic functions. The requirement of absence of divergence u —>■ —00 for 
hnite z leads to the following form for Z^_^i(z): 


f d Vcosh(z-I-C) 


Zn+i{z) = 


C'z"+ 


zdz / z 

1 ( d \ sinh(z -1- Q 

zdz 


n = 2p, 

n = 2p + 1; p = 0, 1...; 


( 9 ) 


1 2 

C = - In —, C = const. 
2 71 
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At first let us to consider the solutions obtained for e and as functions of We 
begin with n = 0, or = -. It is interesting to present the expressions for and ui: 


VI (z) 

e 2 = 


K 


^^^^\Sech\z + C) 
ui{z) = z^(? sech^{z + C), 


( 10 ) 

( 11 ) 


where (10) appears the one-soliton solution exactly the same as the well-known one of the 
Korteweg-de Vries equation. The expressions for n > 1 become more complicated, and one 
can see the solitonic behaviour of and the curves for 'u„^i(z) in Fig.l-Fig.2. 



Fig.l: n = 0, ...6; c = 1; 



Return now to the required variable y. We obtain the explicit form of our solution 

by direct substitution (4) into (9), where — = —c^nin + 1). The resulting formulae are 

c 

complicated and slightly difficult for analytic analysis; it seems to be more convenient to 
present the plots. 

For n = 0 in the function e ^ we have the ’’step” whose altitude depends on the values 
of velocity c and arbitrary constant k. One may see that these curves become higher and 




shift to the right with different rates for rising k. The n i {y) is the positive function whose 
altitude and sharpness of peak depends on c (see Fig.3 - Fig.4). 
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^ 1 {y) 

For n > 1 we can see that the solitonic behaviour of e is retained for different 

values of c and k; the curves become higher, more tight and they shift to the right also with 
increase of c and k. For the cells density the obtained solution has the negative 

section converging to zero for cy —)■ —oo (Figs.5-8). 





The curves for the concentration of the chemical substance are presented in the 

Fig.9. Since n„_,_i(y) have to be positive (nonnegative) we see, that these functions do not 
satisfy this requirement in all domain of dehnition. 



Fig.9: ?;^_^i(y); n = 0,1,...6. Fig.lO: n = 1/5; 8/17; 7; 45. 
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Fig.ll: Uy{y)] v = 7; 45. 


k/C''2=10 

k=l 



In conclusion it seems to be interesting to present the plots for and u^iy) for different 
values of Fig. 10-Fig. 12. It is interesting to see, that there are irregular solutions for 
however the corresponding solutions for Uuiy) are regular (see (6)-(8)). 


III. CONCLUSION 

We investigate the one-dimensional parabolic-parabolic Patlak-Keller-Segel model. One 
of the reductions of this system to ODEs turns out to be integrable. This corresponds to 
the case when the chemical substance diffusion coefficient is equal to 2. After integration 
we obtain the exact soliton solutions in terms of travelling wave variables. The interesting 
fact is that in the plane one of the solution above coincides with the well-known 

Korteweg-de Vries one. However it is very likely that the above solutions are model only and 
they are not useful in practice, and it is interesting to find more useful solutions. Further it 
would be interesting to analyze in more depth the general mathematical properties of this 
system. All these questions require a further investigations. 
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